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Abstract We study the mass spectrum and electromagnetic processes of charmonium system with 
the spin-dependent potentials fully taking into account in the solution of the Schroedinger equation 
I/"") , and the results for the pure scalar and scalar-vector mixing linear confining potentials are compared. It 

is revealed that the scalar- vector mixing confinement is important for reproducing the mass spectrum 
and decay widths and the vector component is found to be around 22%, the long-standing discrepancy 
in Ml radiative transitions of J ftp and ip' is alleviated by means of the state wave functions obtained 
via the Hamiltonian with the full spin-dependent potential. This work also intends to identify few of 
the copious higher charmonium-like states as the cc ones. Particularly, the newly observed X(4160) 
and AT (4350) are assigned as M{2 1 D 2 ) = 4164.9 MeV and M(3 3 P 2 ) = 4352.4 MeV, which strongly 
favor the J PC = 2 h ,2 ++ assignments respectively. The corresponding radiative transitions, leptonic 
and two-photon decay widths have been also calculated for the further experimental study. 
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' 1 Introduction 

O ' 

Due to the impressive increase of experimental results, charmonium (cc) spectroscopy has renewed 
great interest recently, coming along with the striking disagreement with theoretical expectations 
H Hi 9- The unexpected and still- fascinating A (3872) has been joined by more than a dozen other 
charmonium-like states, while the series of vacancy have been left on the cc list. It is urgent to identify 
the possible new members of charmonium family from the abundant observations. 

The QCD inspired potential models have been playing an important role in investigating heavy 
quarkonium, owning to the presence of large nonperturbative effects in this energy region. Most quark 
rS potential models [IS S 0, 1 & El El HI El El EEL Ll6j | have common ingredients under the non- 

relativistic limit, despite some differences in the detailed corrections for relativistic and coupled channel 
effects, which typically are the Coulomb-like term induced by one-gluon exchange plus the long-range 
confining potential expected from nonperturbative QCD. Anyway, the nature of confining mechanism 
has been veiled so far. In the original Cornell model [HEH, it was assumed as Lorentz scalar, which 
gives a vanishing long-range magnetic contribution and agrees with the flux tube picture of quark 
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confinement [iH ] . Another possibility |Qj, 2£i\ 1S that confinement may be a more complicated mixture of 
scalar and timelike vector, while the vector potential is anticonfinning. In pure cc models, the Lorentz 
nature of confinement is tested by the multiplet splitting of orbitally excited charmonium states. 

In addition, numerical precision of the calculation method is important in testing different models. 
As several numerical methods fail in the potentials with 1/r 2 and even higher negative power, the 
0(v 2 /c 2 ) corrections to the quark- antiquark potential has to be usually treated as p erturbation [2ll 
(22|, l20| ■ However, the accuracy of perturbation expansion has been alerted recently [23|, which indicates 
the most significant effect of the different treatments is on the wave functions. The exact solution of 
the full Hamiltonian provides every state with its own wave function, while the perturbative treatment 
leads to the same angular momentum multiplets sharing the identical radial wave function. It is known 
that the radiative transitions, leptonic and double-photon decay widths are quite sensitive to the shape 
of wave function and its behavior at the origin. Each physical particle with different quantum numbers 
should have different state wave function. Thus, it is interesting to study the still-puzzling confining 
mechanism and the numerical precision of the calculation method. 

In our calculation framework, the spin-dependent and -independent interactions have been totally 
taken into account in the Hamiltonian, where the different confining assumptions are compared from 
the mass spectrum and decay properties. The precise wave functions are inspected through the elec- 
tromagnetic processes, i.e. radiative transitions and leptonic decays, which are considered to be a 
nichetargcting test for the overlap of radial integration and the subtle information at the origin. Since 
the relativistic reconstruction of the static confining potential is not unique, it complicates the na- 
ture of confinement. Hence special concerns are focused on the minimal but relatively well- understood 
models, with the aim of gleaning the actual influences of different forms of confinement potential and 
the difference between the exact and the perturbative solutions. 

In the following section, the two potential models are described in detail, along with the adopted 
variational approach and the optimization of parameters. The numerical results are assembled in Secj3l 
a discussion related to the latest experimental results are put in SecHJ Finally, SecE] summaries the 
remarks and conclusions. 



2 Potential models and calculational approach 

The confinement of quarks is assumed to be purely scalar linear type in NR model[2l|, and the scalar- 
vector mixed one in MNR [9(.Once the Lorentz structure of central part are fixed in the two forms, 



NR : V s =br; V v 
MNR : V s =b(l- e)r; V v 



4as 
3 r ' 



4«s 
3 r 



ebr, 



(1) 
(2) 



the spin-orbit term and the tensor term can be directly derived from the standard Breit-Fermi ex- 
pression to order (v 2 /c 2 ) with the charmed quark mass m c . Summarily, the interaction potentials 
are 



Vnr = — r — 
3 r 



V, 



MNR 



h br 

S r 



b r 

32na s 
9m 2 



32ira s ~ 

d a (r)S c ■ S- c 



9m 2 
<5> CT (r)S c ■ S 



2a s 



' 2a s b 
m 2 r 3 2m 2 r 

{As- 1)6 



L S 



2m 2 r 



L S 



sb 



3m 2 r 3 I2m 2 7 



(3) 



T, (4) 



where L is the orbital momentum and S is the spin of charmonium. In the mixed- confining model, 
e stands for the vector exchange scale. The singularity of contact hyperfine interaction within the 
spin-spin term has been smeared by Gaussian as in Ref.iH, 6„{r) = O/0F) 3 &-^ r \ The involved 
operators are diagonal in a | J, L, S) basis with the matrix elements, 



(S c -S e ) = is(S+l)- j, 
(L-S) = \[J{J+l) -L(L+1)- S(S+1)}, 



(5) 
(6) 
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Table 1 Basis-related parameters for the potential models 





Parameter 


NR 


MNR 






10 


8 


Basis Space 


n [GeV" 1 ] 


0.4 


0.4 




rn m „ [GeV" 1 ] 


15.8 


6.8 




m c [GeV] 


1.4786 


1.5216 






0.5761 


0.6344 


Potential Model 


b [GeV 2 ] 


0.1468 


0.1361 




a [GeV] 


1.1384 


1.2058 




e 




-0.2193 



r(S c T)(S g -r)-(S c -S g ) 



6 ((L • S)) 2 + 3(L ■ S) - 25(5 + 1)L(L + 1) 



(7) 



6(2Z- l)(2L + 3) 

Instead of separating the spin-dependent interactions into leading order portions, we solved the Schroedinger 
equation of the unperturbed Hamiltonian with full potentials, 



r— V 2 + K e (r) 



E 



4>(t) = 0, 



(8) 



where fiR — m c /2 is the reduced mass. Here, V c£ (r) is Vnr or Vmnr which includes the spin- 
independent interactions as well as the spin-dependent ones. The Hamiltonian includes the full po- 
tential, enable us to maintain the subtleness of wave function. With the help of a well-chosen set of 
Gaussian basis functions, namely Gaussian Expansion Method [24| , the singular behavior of 1/7- 3 in 
spin-dependent terms at short distance can be refined variationally. Then, the wave functions ipim(r) 
are expanded in terms of a set of Gaussian basis functions as 



Tlmax 

E 

n=l 



C n l 



2U„ 



^") 2 y ;m (f), 



(9) 



(10) 



The basis-related parameters {n max ,ri,r n } are determined by the variational principle, resulting 
the reasonably stable eigen solutions and differ slightly according to the explicit forms of potential 
models. Table Q] shows the utilized basis for the two potential models. 

To determine the parameters (a = a%, 0,2, • • • , ap) appearing in the potentials, a merit function \ 2 
has been defined to search the best-fit parameters by its minimization, 



X 



N 

E 

i=l 



Mf 



(11) 



where N denotes the number of targeted data, the <r are the associated errors and M, represents 
the experimental and theoretical values. Given a trial set of model-depended parameters, a procedure 
calculating the \ 2 ( a ) are developed to improve the trial solution with the increments <5a and repeated 
until x 2 (a+(5a) effectively stops decreasing. The increments <5a are solved by the set of linear equations 



with 



2 da k 



N 1 

ak i = E 



i=l 



dMf h (a) dM} h (a 



dai 



9a, 



(12) 

(13) 
(14) 
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The iteration has been taken in the numerical analysis of nonlinear systems. Here, the partial derivative 
has to be solved numerically because the M J " l (a) is not analytic. 



3 Numerical Results 

3.1 Mass Spectrum 

Combining the nonrclativistic kinetic term and the interactions terms, we diagonalize the full Hamilto- 
nian, which leads to a generalized eigenvalue problem. After the fitting procedure, the energy levels and 
corresponding eign functions of the two potential models are obtained with its optimized parameter 
set (shown in Table [1]). 

The vector component of linear confining are fitted to be 21.93% for the MNR model. This result 
consists with the one- fifth vector exchange given by Ref. 23] . If defined an additional r-independent 
adjustable parameter in the vector linear confinement, the value of mixing coefficient should be — 1 
for some decay considerations [20]. In addition, the scalar- vector mixing model requires less linear 
potential slope, but more splitting scale from spin-spin interaction. 

Totally 60 states have been calculated with the two potential models, and summarized in Table [5J 
We fit the mass of well-established states marked by *. Except somewhat low at Xo(l 3 Po), both of the 
two potential models are overall good to reproduce the spectrum. 



3.2 Leptonic Decays 

The lowest-order expressions of electronic decay width the first-order QCD corrections [26| are 

4a 2 e 2 2 / 16a s (m c )\ 

r ee(nS) = | R nS {0) | 2 \l - y ~ J , (15) 

25« 2 e 2 2 / 16 a s (m c ) \ 

F ^ nD) =2M^M* lR " D{0)l ^-Y—^J' (16) 

where e c is the c-quark charge in units of | e |, a = 1/137.036 is the fine-structure constant, M n s and 
M n D are the (n r + l)th S-wave and D-wave state mass respectively with the radial excitation number 
n r . Note that here a s (m c ) and the a s are essentially the strong coupling constants of different mass 
scales, and we adopt a s (m c ) = 0.26 as in Ref. [13, [28j]. R n s{0) is the radial 5* wave function at the 
origin, and R'^ D (0) is the second derivative of the radial D-wave function at the origin. Within the 
Gaussian basis space, the analytic formula for R(0) is explicitly presented in the APPENDIX. Table |3] 
compares the leptonic decay given in Ref. [13] with our results in NR and MNR models. 



3.3 Two-photon Decay 



The two-photon decay widths are important to identify the potential charmonium states, i.e. A(3915), 
Z(3930) etc. and the latest A (4160), A (4350). With the first-order QCD radiative corrections [2fJ, the 
two-photon decay widths of 1 S'o, 3 Po and 3 i"2 explicitly are 



p 77 ( 5o) 



r 77 ( 3 p ) — 



3a 2 e4 | R nS (0) | 2 
m 2 

27a 2 e^ | R' nP (0) | 2 



r 77 ( 3 p 2 ) 



ml 
36a 2 ef 



a s (m c ) f 7T 2 
7T ^ 3 

a s {m c ) ( 7T 2 
3 



20 
3 

28 
" ~9 



5ml 



1 



16 a s (m c ) 



(17) 
(18) 
(19) 



where R' nP (0) is the first derivative of the radial P-wave function at the origin. Our numerical results 
of two-photon decay widths are shown in Table 0] in comparison with Refs. 28, 30]. 
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Table 2 Experimental and theoretical Charmonium mass spectrum. The masses are in units of MeV, and 
the * denotes the states used in the optimization of potential parameters. Our full-potential calculation results 
are listed in comparison with the perturbative results of NR model and its relativized extension GI model [2j . 



State Expt. [25] Ref. [21] Ours 

NR GI NR MNR 



J/^(l 3 5i) 
V'(3 3 Si)* 

^(4? Si)* 
X2(1 3 P 2 )* 

xo(i 3 p y 

hci^PxT 

X2(2 3 P 2 ) 

Xi(2 3 Pi) 

Xo(2 3 Po) 

ft c (2 1 i\) 

X2(3 3 P 2 ) 

Xi(3 3 Pi) 

Xo(3 3 Po) 

hep 1 ft) 

X2(4 3 P 2 ) 
Xi(4 3 Pi) 
Xo(4 3 Po) 
M^Pi) 

^(l 3 Pi)* 

M2 3 D 2 ) 
^(2 3 J Di)* 
Vc2(2 1 D 2 ) 

M3 3 D 2 ) 

<K3 3 L>i) 

V a 3(3 1 D 2 ) 

X4(1 3 P 4 ) 
X3(1 3 P 3 ) 
X2(1 3 P 2 ) 

M^s) 

X4(2 3 P 4 ) 

X3(2 3 P 3 ) 

X2(2 3 P 2 ) 

hc3(2 1 F 3 ) 

X4(3 3 P 4 ) 

X3(3 3 P 3 ) 

X2(3 3 P 2 ) 

hc3(3 x F 3 ) 

^5(1 3 G 5 ) 
^ 4 (1 3 G 4 ) 
^ 3 (1 3 G 3 ) 



2980.3 ± 1.2 


2982 


2975 


2990.4 


2978.4 


3096 916 ± Oil 


3090 


3098 


3085.1 


3087.7 


3637 ± 4 


3630 


3623 


3646.5 


3646.9 


3686 09 ± 04 


3672 


3676 


3682.1 


3684.7 




4043 


4064 


4071.9 


4058.0 


4039 ± 1 


4072 


4100 


4100.2 


4087.0 




4384 


4425 


4420.9 


4391.4 


4421 ± 4 


4406 


4450 


4439.4 


4411.4 


occfi on 4- n no 


oOOD 


OOOU 


O001.4 


oOOy.o 


3510 66 ± 07 


3505 


3510 


3500.4 


3517.7 


3414 75 + 31 


3424 


3445 


3351.9 


3366 3 


3525.42 ± 0.29 


3516 


3517 


3514.6 


3526 9 




3972 


3979 


3979.8 


3973.1 




3925 


3953 


3933 5 


3935 




3852 


3916 


3835.7 


3842.7 




3934 


3956 


3944.6 


3941.9 




4317 


4337 


4383.4 


4352.4 




4271 


4371 


4317.9 


4298.7 




4202 


4292 


4216.7 


4207.6 




4279 


4318 


4333.9 


4309.7 








4736.7 


4703.1 








4620.3 


4590.5 








4551.8 


4521.7 








4639.5 


4606.7 




3806 


3849 


3814.6 


3812.6 




3800 


3838 


3807.7 


3820.1 


3772 92 ± 35 


3785 


3819 


3785.3 


3808 8 




3799 


3837 


3807.3 


3815.1 




4167 


4217 


4182.9 


4166.1 




4158 


4208 


4173.7 


4168.7 


4153 ± 3 


4142 


4194 


4150.4 


4154.4 




4158 


4208 


4173.7 


4164.9 








4572.5 


4526.5 








4558.8 


4523.6 








4525.8 


4502.2 








4559.7 


4521.4 




4021 




4037 4 


4024.7 




4029 


4097 


4044.0 


4047.6 




4029 


4092 


4042.4 


4059.7 




4026 


4094 


4041.1 


4040.8 




4348 


4425 


4371.1 


4344.7 




4352 


4426 


4374.4 


4362.4 




4351 


4422 


4369.9 


4369.8 




4350 


4424 


4372.3 


4356.8 








4744.4 


4684.8 








4747.7 


4698.5 








4743.9 


4704.2 








4745.9 


4694.3 




4214 


4312 


4236.9 


4213.5 




4228 


4320 


4250.6 


4244.7 




4237 


4323 


4258.2 


4267.4 




4225 


4317 


4247.1 


4237.9 
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Table 3 Leptonic decay width, in units of keV. 



Particle State 



r; 



Ref. [27] 



Ref. [29] 



NR 



MNR 



r; 







r. 



Expt.[25j 



■0(4040) 
0(4415) 
V>(3770) 
■0(4160) 



2 3 Si 

3 3 ,Sl 

4 3 Si 
l 3 Di 
2 3 L»i 



11.8 
4.29 
2.53 
1.25 
0.055 
0.066 



6.6 
2.4 
1.42 
0.7 
0.031 
0.037 



12.13 
5.03 
3.48 
2.63 
0.056 
0.096 



5.6 
2.3 
1.9 
1.3 
0.089 
0.15 



3.1 
1.3 
1.0 
0.70 
0.050 
0.084 



6.0 
2.2 
1.8 
1.3 
0.079 
0.13 



3.3 
1.2 
0.98 
0.70 
0.044 
0.073 



5.6±0.14±0.02 
2.4±0.07 
0.86±0.07 
0.58±0.07 
0.27±0.018 
0.83±0.07 



Table 4 Two-photon decay width, in units of keV. 



State 


Ref. [28] 


Ref. [30] 


NR 


MNR 


Expt. [25] 


l'So 


5.5 


10.94 


7.4 


7.5 


6 7 +0 ' 9 

D -' -0.8 


2 1 5o 


1.8 




3.2 


2.9 








2.9 


2.5 




4 1 S'o 






2.0 


1.8 




l 3 Po 


2.9 


6.38 


11 


10.8 


2.29±0.18 


2 3 P 


1.9 




7.7 


6.7 




3 3 P 






7.9 


6.5 




1 3 P 2 


0.50 


0.57 


0.29 


0.27 


0.50±0.03 


2 3 P 2 


0.52 




0.43 


0.39 




3 3 P 2 






0.81 


0.66 





3.4 Radiative Transition 

Because radiative transition is sensitively dependent on the detailed features of the wave functions, it 
is of great interest to have careful inspection. Besides, it has been noticed that the known Ml rates 
showing serious disagreement between the previous theoretical calculation [2lT | and experiment. In our 
full-potential calculation framework, the wave functions are directly associated with the eigenvectors of 
Hamiltonian corresponding to the masses. The El transition rate between an initial charmonium state 
i of radial quantum number nj, orbital angular momentum Li, spin Si, and total angular momentum 
Ji, and a final state / is given by Ref. [3l[ as 

r E1 (nf' +1 L, j4 H- n f Sf+1 L fjf ) = -C fi 6 SiSf e 2 c a\(i> f \ r \^)\ 2 E 3 y -^, (20) 



In the above formulas, e c is the charge of c-quark in unit of | e |, and Mi,Ef represent the eigen mass of 
initial state and the total energy of final state respectively. The momentum of the final photon equals 
E y = (Alf — Mj)/(2Mi) in the nonrelativistic approximation (32|. The angular matrix element C/j is 

C fi = max(L t , L f )(2J f + 1) | L J_ J \ } . (22) 

The Gaussian expanded wave functions give rise to the analytic formulas for the overlap integral, 
i.e. Eq. (|A.8j) . and the transition matrix elements Eq. (|A.10|) . The n max stands for the dimension of 
Gaussian basis as defined above. The numerical results of El transition rates are presented in Table [S] 
and Table HI as well as those of Ml transition in Table [7] 
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Table 5 2S, 3S, IP and 2P El radiative transitions 









Ref.[21] 






Ours 




P Tori ru^T/l 


T 1 

Initictl H16SOI1 


r inal meson 




McV] 


Pthy 


[Ke v j 




[MeV] 


Pthy 


[Ke v j 








NR 


CtI 


INK, 


GI 


NR 


TV /TT\TT> 


INK 


MNR 




^'(2 3 S 1 ) 


X2(1 3 P 2 ) 


128 


128 


38 


24 


128 


123 


43 


39 


26.6 ± 1.1 




Xi(l 3 Pi) 


171 


171 


54 


29 


177 


163 


48 


38 


28.0 ± 1.2 




Xo(l 3 Po) 
^(l'Pi) 


261 


261 


63 


26 


315 


305 


34 


29 


29.2 ±0.9 




111 


119 


49 


36 


130 


118 


72 


56 




X2(1 3 P 2 ) 


455 


508 


0.70 


12.7 


512 


494 


6 


6.44 


< 1.4 




Xi(l 3 Pi) 


494 


547 


0.53 


0.85 


556 


530 


0.76 


0.44 


< 0.9 




Xo(l 3 Po) 


577 


628 


0.27 


0.63 


680 


657 


8.9 


8.1 






485 


511 


9.1 


28 


519 


496 


6.1 


7.7 




X2(1 3 P 2 ) 


J/i>{i z Si) 


429 


429 


424 


313 


436 


440 


421 


405 


384.2 ± 16 


Xi (i 3 Pi) 




390 


389 


314 


239 


391 


404 


330 


341 


295.8 ± 13 


Xo(l 3 Po) 




303 


303 


152 


114 


256 


267 


97 


104 


119.5 ±8 


hci^Pi) 


^'(2 3 Si) 


504 


496 


498 


352 


485 


506 


465 


473 




X2(2 3 P 2 ) 


276 


282 


304 


207 


287 


278 


300 


264 




Xi(2 3 Pi) 




232 


258 


183 


183 


243 


242 


243 


234 




Xo(2 3 Po) 
h^Pi) 




162 


223 


64 


135 


151 


155 


77 


83 




r i i /~\ 3 a \ 

j/ip{i s 


285 


305 


280 


218 


287 


284 


297 


274 




X2(2^P 2 ) 


779 


784 


81 


53 


794 


787 


108 


111 




Xi(2 3 Pi) 




741 


763 


71 


14 


757 


756 


27 


33 




Xo(2 3 Po) 
h^Pi) 




681 


733 


56 


1.3 


677 


681 


30 


28 






839 


856 


140 


85 


839 


846 


104 


116 




X2(2 3 P 2 ) 




163 


128 


88 


29 


162 


157 


81 


76 






2 (1 3 £> 2 ) 


168 


139 


17 


5.6 


168 


150 


14 


10 






V(i 3 £>i) 


197 


204 


1.9 


1.0 


190 


161 


0.98 


0.64 




Xi(2 3 Pi) 




123 


113 


35 


18 


124 


113 


37 


30 






0(l 3 £>l) 


152 


179 


22 


21 


145 


124 


16 


11 




Xo(2 3 Po) 

frc^Pl) 


0(l 3 £>l) 


81 


143 


13 


51 


50 


34 


4.2 


1.4 




mc(i 1 D 2 ) 


133 


117 


60 


27 


135 


125 


61 


51 





4 Discussion 



4.1 V(3686), 0>(377O), 0(4160) 



The world-average r ee (ip(3770)) is 0.265 ± 0.018 keV 25]. The significant leptonic width implies that 
there is a sizeable S-D mixing between the 2 3 Si and l 3 £>i state, since it is expected to be highly 
suppressed if ^(3770) is a pure D-wave state (shown in Table The mixing arises both from the usual 
relativistic correction terms and coupling to strong decay channels, and will affect the El transition 
rates of ■0(3770) and ■0(3686) (32j . We investigate the mixing under the present calculation framework 
from two scenarios respectively: electronic annihilation and dipole transition. Assuming the ■0(3686) 
and ?/>(3770) to be a mixture of a l 3 .Di and a 2 3 5i state as 



Then, Eqs. ([TBI I16f can be expressed as 



V>(3686)) = | 2 3 Si)cos6>- | l 3 £>i)sin0, 
0(3770)) = | 2 3 Si)sin(9+ | l 3 L>i)cos6». 



(23) 



r ee (V>(3686)) 



r ee (0(377O)) 



Ml, 
Aa 2 e 2 c 



R 2S (0)cos9 



R 2S (0) sm6- 



2y/2m 2 c 



R'{ D (0)ahx9 



: R'{ D (0) COS ( 



(24) 



(25) 



Through the fitting of the ratio P"(-0(377O) — * e+e )/-P"(0(3686) — > e + e ), combining the exper- 
imental data j25| and the radial wave function values at the origin, i.e. R 2S (0) = 0.556 GeV 3 / 2 and 
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Table 6 3P, ID and 2D El radiative transitions 









Ref. [21] 






Ours 




1 expt. 




nidi 111L.0011 


JJJ<y 


[MeV] 


r t h y 


|Kc V J 




[MeV] 


Pthy 


[Kc V J 


I V J 






NR 




INK 


GI 


NR 


MIN K 




MNR 




X2(3 3 P 2 ) 


^(3 3 5i) 


268 


231 


509 


199 


273 


257 


434 


360 




Xi(3 3 Pi) 




225 


212 


303 


181 


212 


206 


292 


272 




Xo(3 3 Po) 




159 


188 


109 


145 


115 


119 


63 


72 




V/(2 3 Sr) 


229 


246 


276 


208 


254 


244 


377 


332 




X2(3 3 P 2 ) 


585 


602 


55 


30 


645 


617 


112 


97 




Xi(3 3 Pi) 




545 


585 


45 


8.9 


589 


570 


30 


31 




Xo(3 3 Po) 




484 


563 


32 


0.045 


501 


490 


25 


19 




h^Pi) 




593 


627 


75 


43 


633 


612 


95 


90 




X2(3 3 P 2 ) 


J/V(l 3 50 


1048 


1063 


34 


19 


1106 


1081 


62 


59 




Xi(3 3 Pi) 




1013 


1048 


31 


2.2 


1057 


1040 


10 


11 




Xo(3 3 P ) 
MS'Pi) 




960 


1029 


27 


1.5 


980 


971 


23 


22 






1103 


1131 


72 


38 


1135 


1126 


73 


76 




X2(3 3 P 2 ) 


</>3(2 3 P> 3 ) 


147 


118 


148 


51 


196 


182 


279 


231 






</> 2 (2 3 P> 2 ) 


156 


128 


31 


9.9 


205 


180 


48 


34 






^(2 3 P>!) 


155 


141 


2.1 


0.77 


227 


194 


3.2 


2.1 




Xi(3 3 Pi) 


</> 2 (2 3 P> 2 ) 


112 


108 


58 


35 


142 


128 


121 


92 






^(2 3 Di) 


111 


121 


19 


15 


164 


142 


50 


34 




Xo(3 3 P ) 


V(2 3 £>i) 


43 


97 


4.4 


35 


66 


53 


21 


11 






V2c(2 1 D 2 ) 


119 


109 


99 


48 


157 


142 


205 


158 




X2(3 3 P 2 ) 




481 


461 


0.049 


6.8 


532 


506 


1.6 


1.4 






^ 2 (1 3 -D 2 ) 


486 


470 


0.0091 


0.13 


538 


500 


0.76 


0.42 






</>(l 3 P>i) 


512 


530 


0.00071 


0.001 


557 


510 


0.11 


0.062 




Xi(3 3 Pi) 


</> 2 (l 3 A0 


445 


452 


0.035 


4.6 


480 


452 


1.2 


1.8 






V>(l 3 #i) 


472 


512 


0.014 


0.39 


500 


462 


0.11 


0.00029 




Xo(3 3 Po) 


</>(l 3 I>i) 


410 


490 


0.037 


9.7 


409 


380 


28 


30 








453 


454 


0.16 


5.7 


495 


466 


0.21 


0.56 






X2(1 3 P 2 ) 


242 


282 


272 


296 


254 


245 


340 


302 






X2(1 3 P 2 ) 


236 


272 


64 


66 


248 


252 


79 


82 






xi(l 3 A) 


278 


314 


307 


268 


295 


290 


321 


301 




w(l D\) 


X2(l P 2 ) 


208 


208 


4.9 


3.3 


227 


241 


6.8 


8.1 


< 24.6 




Xi(l 3 Pi) 


250 


251 


125 


77 


274 


280 


146 


153 


79.2 ± 16 




Xo(l 3 Po) 


338 


338 


403 


213 


409 


417 


367 


362 


199.3 ±25 




^(l'Pi) 


264 


307 


339 


344 


281 


277 


398 


374 






X2(1 3 P 2 ) 


566 


609 


29 


16 


584 


563 


23 


23 






X2(1 3 P 2 ) 


558 


602 


7.1 


0.62 


576 


565 


0.85 


1.5 






xi(l 3 A) 


597 


640 


26 


23 


619 


600 


38 


41 




V>(2 3 P>i) 


X2(1 3 P 2 ) 


559 


590 


0.79 


0.027 


556 


552 


0.18 


0.080 


< 1.3 




Xi(l 3 A) 


598 


628 


14 


3.4 


599 


588 


3.6 


5.5 


< 0.7 




Xo(l 3 Po) 


677 


707 


27 


35 


722 


713 


99 


109 




h c {2 1 D 2 ) 


^(l'Pi) 


585 


634 


40 


25 


607 


589 


38 


41 





R'{ D (0) = 0.144 GeV 7 / 2 , we find two solutions for the mixing angle: 

(9 = 8° ±0.4° or 9 = -30° ±0.4°. (26) 

These solutions are smaller than (12±2)°or— (27±2)° determined by leptonic decay widths in (35l.[27j. 
partially due to more precise measurements. Impending, we illustrate the mixing-dependency of El 
transitions of our model in Fig. [T]and Fig. [2] The decay widths of ^(3686) —> xxco, 7Xci obviously 
favor 8 — 8°, while the matching angle for 7Xc2 locates around 9 = —30° and 9 = 50°. The 9 = 50° 
would be too large for the two states mixing. On the other hand, a positive and slight mixing is 
compatible for ?/>(3770) — > 7Xc0, 7Xc2 consisting with 9 — 8°, except the challenging -fXa fitting with 
a small but negative angle. This discrepancy reveals that S-D wave mixing would be important for 
i/>(3770) and ^(3686), but not sufficient. The mixing-angle functions for the two out-standing decays 
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Table 7 Ml radiative partial widths 



Initial 


r mal 




Ref. [211 








Durs 




1 expt. [KeV J 


niGSOii 


niGson 




MeV] 


Pthy 


[Ke V J 




IJYLc V 1 


1 thy 


|Ke v j 






NR 




IN rt 


GI 


NR 


MNR 


NR 


MNR 




J/M1 3 S X ) 
0'(2 3 Si) 


^(l'So) 
V^So) 


116 


115 


2.9 


2.4 


93 


107 


1.5 


2.2 


1.58 ±0.37 [25] 


48 


48 


0.21 


0.17 


35 


38 


0.086 


0.096 


0.143 ±0.027 ±0.092 [33] 




639 


638 


4.6 


9.6 


627 


639 


3.1 


3.8 


0.97 ±0.14 [34] 




J/0(l 3 Si) 


501 


501 


7.9 


5.6 


518 


516 


6.1 


6.9 


^(3 3 Si) 


29 


35 


0.046 


0.067 


28 


29 


0.043 


0.044 






382 


436 


0.61 


2.6 


429 


416 


0.70 


0.71 






0'(2 3 Si) 


922 


967 


3.5 


9.0 


960 


958 


3.2 


3.7 






312 


361 


1.3 


0.84 


371 


356 


1.7 


1.6 






J/V(1 3 S!) 


810 


856 


6.3 


6.9 


867 


854 


5.9 


6.5 






X2(1 3 P 2 ) 


360 


380 


0.071 


0.11 


374 


364 


1.3 


1.2 






Xi(l 3 Pi) 


400 


420 


0.058 


0.36 


419 


401 


0.16 


0.13 






Xo(l 3 Po) 
h e {l 1 I\) 
^(l'Pi) 


485 


504 


0.033 


1.5 


548 


534 


5.6 


5.3 




X2(2 3 P 2 ) 


430 


435 


0.067 


1.3 


438 


421 


1.0 


0.89 




Xi(2 3 Pi) 


388 


412 


0.050 


0.045 


397 


387 


0.15 


0.13 




Xo(2 3 Po) 




321 


379 


0.029 


0.50 


308 


303 


4.8 


4.8 




X2(2 3 P 2 ) 
Xi(2 3 Pi) 
Xo(2 3 Po) 
X2(1 3 P 2 ) 
Xi(l 3 Pi) 










340 
382 
470 
712 
753 


323 
358 
442 
685 
719 


1.7 
0.17 
3.5 
2.3 
0.27 


1.5 
0.12 
2.7 
2.1 
0.22 






Xo(l 3 Po) 
hc&Pi) 

h c (2 1 P l ) 

hri^Pl) 










871 


840 


8.9 


8.2 




X2(3 3 P 2 ) 










417 


391 


1.2 


0.99 














783 


747 


1.8 


1.5 




Xi(3 3 Pi) 










357 
729 


342 
703 


0.16 
0.22 


0.12 
0.20 




Xo(3 3 Po) 


h c {2 1 Pi) 
^(l'Pi) 










263 


257 


3.7 


3.3 














644 


626 


5.3 


5.1 





may be modified by the first-order relativistic corrections to their wave functions, since the relativistic 
corrections play an essential role in El radiative transitions indeed, as suggested in [35j . 

The S-D mixing may be more serious for the higher charmonium above the open-charm thresh- 
old [H, [H, [13, H3- As shown in Table El the experimental leptonic decay width of the -0(4160) is 
significantly larger than the potential model predictions as 2 3 Di. In Ref. [33|, it is explored that the 
possible mixing in (n + l) 3 Si-n 3 Di determined by the ratio of di-electron widths and get a large 
mixing angel ~ 34° (n = 2,3), while the hyperfine effects are neglected in the wave functions by a 
universal potential. In the current framework, the spin-dependent interactions at 0(v 2 /c 2 ) are thor- 
oughly considered to determination of the state wave functions. Besides the distinct leptonic decay 
widths, it is found that the theoretical El radiative rates of the pure 2 3 Z?i and 3 3 Si are incomparable 
with the available experimental data. This discrepancy seems cannot be eliminated by the relativistic 
corrections, hence dissimilar with the situations of 0(3770) and ^(3686). We extend above pattern to 
the 0(4O4O)- , 0(416O) and ip(4160)-ip(4Al5) mixing scenarios, and also independently via the leptonic 
decays as well as their El radiative rates. Combining the corresponding values of wave functions at 
the origin, i.e. R 3S (0) = 0.556 GeV 3 / 2 , R 4S (0) = 0.508 GeV 3 / 2 and R% D (0) = 0.201 GeV 7 / 2 and the 
central values of experimental data (4040)) = 0.86 ± 0.07 keV, f ee (^(4160)) = 0.83 ± 0.07 keV 
and Pee (0(4415)) = 0.58 ± 0.07 keV [25|], the mixing angles determined by the leptonic decay widths 
are listed in Table [5] 

In FigO one can figure out that the mixing angles given by the P ee , i.e. 29.8° or —60.7°, are 
still acceptable for the El transitions. Then, it is in a dilemma that the P(0(4O4O) — > JXd) and 
P(0(416O) ->■ 7x c2 ) favor 9 = 29.8°, while r(^>(4040) 7x c2 ) and P(0(416O) -> 7Xci) support 9 = 
—60.7°. Both of the two solutions seems to be equally important, and it is difficult to rule out any one of 
them from the current perspective, possibly due to the mixing model being oversimplified. Nevertheless, 
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Fig. 1 The El radiative transition rate of t/j(3686) — >• 7l 3 Pj dependent on the various mixing angle 9. Because 
of the relatively precise and very close measurements for the three channels, we illustrate the central values of 
r exp t. 25] with single lines rather than the error-associated areas. 



350 




-90° -60° -30° 0° 30° 60° 90° 



Fig. 2 The El ra diat ive transition rate of ^(3770) — > jl 3 Pj dependent on the various mixing angle 9. The 
experimental data [25[ are illustrated by the central values as well as the related errors. 

Table 8 The 3 3 Si-2 3 Di and 4 3 Si-2 3 Di mixing for the higher charmonium di-electron decays. Our predictions 
are obtained with the MNR potential model. Since the inverse definitions about mixing angle, these assignments 
in Ref. [3g] are reasonably exchanged here for a proper comparison. 



Mixing states 




Mixing angle 






Ref. [36] 


Ref. [37] 


Ours 


3S-2D 


357-55° 


34.8°/ - 55.7° 


29.8°/ -60.7° 


4S-2D 


24.8°/ - 58.5° 
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Fig. 3 The mixing angle ^-dependent El radiative transition rates of ^(4040) — > 7Xci (solid), ^(4040) — > 7Xc2 
(dash), t/j(4160) — > 7Xci (short dash dot) and i/>(4160) — > 7Xc2 (short dash). The experimental upper limits 
|25lj are represented in the same line-shape with the corresponding channels. 



it would be still interesting to glean something for the further understandings and revisiting with more 
precise measurements. 



4.2 A(3872) and X(3915) 

The first surprising charmonium-like state A (3872) [381 ] is still far away from sufficiently understanding. 
Its quantum numbers have not been decided yet from J PC = 1 ++ and J PC = 2 h [39[. The radiative 
decay jJ/ip [IS EH bolsters the 1 ++ assignment, while the J p — 2~ is ratherish favored by a study of 
the 3ir invariant mass distribution in J/ipuj decays [42l l2j|. Assuming the charmonium structure, these 
J PC imply the possible candidates, i.e. Xci(2 3 -Pi) and ^(l 1 -^). As the spectrum listed in Tabled! the 
masses of the two states remain a gap around 50 MeV with the world-average mass M = 3871.56±0.22 
MeV [25]. Combining the observed branching fraction |40j and the full width, it is evaluated that the 
radiative decay widths are r(A(3872) -> J/i(jj) > 20.7 keV and r(A(3872) -> ip(2S)i) > 0.069 
MeV. The theoretical results of the radiative decays for the conventional charmonium candidates are 
apparently too large, and the radio of the two rates is nearly two times larger than the observation. 
Hence, the cc nature of A (3872) seems unlikely in the present framework. 

Another interesting resonance named A(3915) has been announced by Belle Collaboration (43[, 
which is found in the processes 77 — > J/ipw. A(3915) clearly has the positive C-parity, but the J p 
assignment remains to be determined. Owning to the same production and decay signatures, the 
V(3940) pi. |45| and the A(3915) are considered to be the same state prevailingly, which is close to 
those of the Z(3930). If the Z(3930) (discussed in the coming section) occupies the assignment of x' C 2 
@,|46|, the possibility of A (3915) can be ruled out, or else they are the same state as Xc2(2 3 P2). The 



world-average mass and width of the A(3915) are M = 3917.4 ± 2.7 MeV, F = 28±^° MeV |25|. In 
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our predictions, the mass of Xc2(2 P2) is about 60 MeV higher than the experimental data for the 
A(3915). Even if the mass is roughly plausible, the conventional cc structure of A(3915) would likely 
encounter trouble in accommodating the product of the two-photon decay width and the branching 
ratio to J/ipiv, i.e. T( 77 ) x Z?(A(3915) -> J/ipu) = 61 ± 17 ± 8 eV for J p = 0+ and 18 ± 5 ± 2 
eV for J p = 2+. As shown in Table HI the two-photon decay width is T 77 (2 3 P 2 ) = 0.39 keV in the 
MNR potential model. Based on the central value of the world-average width, P(A(3915) — > J/i/jcu) is 
expected to be 1.29 MeV ( J PC = 2+) or 4.38 MeV (J PC = 0+) approximately, which is over ten times 
larger than that of the typical charmonium. On the other hand, the charmonium candidate Xco(2 3 ^o) 
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with J p = ++ may be slightly better in 1^(77), but worse in mass spectrum. The strong-decay 
inspection of the X(3915) also disfavor the structure as Xco(2 3 Po) [47j, [48|. It has been interpreted as 
the multiquark and molecule state [49|, d(| HJ- So far, the nature of JT(3915) is still murky, and the 
charmonium system have seemed inefficient to provide a consensus of interpretation. 

4.3 V(3940) and Z(3930) 

The X(3940) was seen by the Belle Collaboration in the recoiling spectrum of J/tp in the e + e~ 
annihilation process [H2|. The mass and width of V(3940) are M = 3942±£ ± 6 MeV, P = Z7±f 5 ± 8 
MeV [25]. This state is not seen to decay into DD but does decay into DD* , suggesting that it has 
unnatural parity and can be a candidate for the rjcffiSo). The predicted mass is much higher by the 
potential models (2^, H3, HO, HH , including the present work. Taking the coupling channel effect into 
account [H[, the induced spin-spin splitting between 3 3 5i and 3 1 5o will be increased and consequently 
lower down the mass of state. We predicate the radiative transitions and two-photon decay width 
of 3 1 So for testing the nature of V(3940) as the charmonium candidate. 

The Belle Collaboration has observed Z(3930) state [5J| in two-photon fusion, 77 — > Z(3930) — > 
DD. The world-average mass and full width are M = 3927.2 ± 2.6 MeV, P = 24 ± 6 MeV [25|. And 
the theoretical prediction of mass within our framework is almost the same as the perturbative results, 
nearly 40 MeV higher than the observation. Due largely to the agreement with the detected angular 
distribution [5J, [55| for J = 2 and to the plausible mass of a 2P cc state, Z(3930) is now widely 
accepted as Xc2(3930) [H, [13, [H, 0, H . The value of the partial width P( 77 ) x £(£(3930) -> DD) = 
0.24 ± 0.05 ± 0.04 keV was found in the Babar experiment [HI]. Our prediction of P( 77 ) is 0.39 keV 
and 0.43 keV within the MNR model and the NR model respectively. With the theoretical open-charm 
strong decay width P(Z(3930) -> DD) = 21.5 MeV [El], one can directly evaluate P( 77 ) « 0.27±0.07 
keV. Our prediction of the two-photon width for Xc2(2 3 P2) is reasonably consistent with this estimation, 
supporting Z(3930) assigned as charmonium 2 3 p2. 

4.4 X(4160) 

The Belle Collaboration has reported the V(4160) in the process e + e~ — > J/ipD*D* with a significance 
of 5.1 a [58]. Its mass and width are M = 4156^^ ± 15 MeV and r = 139^^^ ± 21 MeV respectively. 
Since within this process e + e~ annihilate into one photon, the charge parity of V(4160) should be 
C = +. The possible interpretation of its structure from the view of the production rate in e + e~ — > 
J/V>V(4160) favors the candidate Y r n(3 3 Pn) [Hj], although the predicted mass differs around 30 MeV in 
the color screened potential model [27]. In our calculated spectrum shown in Tabled our prediction of 
r) c {2 1 D 2 ) is M(2 1 D 2 ) = 4164.9 MeV in MNR model and 4173.7 MeV in NR model. It is suggested that 
V(4160) can be excellently reproduced by the charmonium spectrum as 77 C (2 1 Z?2) in our calculation 
framework. Investigating the strong decay properties of V(4160), Ref. [47| has indicated that the decay 
width of i] c (2 1 D2) agrees well with the Belle data of V(4160), which includes the absence of theoretical 
forbidden channel 77 c (2 1 £> 2 ) — > DD. Other possible candidates, i.e. Xco(3 3 Po), Xci(3 3 Pi),^ c (4 1 5'o) are 
disfavored from this perspective. The charmonium structure of V(4160) as r/ c (2 1 Z?2) can be further 
explored by electromagnetic transition rates, which have been predicted in this paper as well. 

4.5 V(4350) 

The V(4350) found in the invariant mass spectrum of J/ip(f> by the Belle Collaboration [l^], is of 
M = 4350^5 j ± 0.7 MeV and T = ISStlY ± 4.1 MeV. Its assignment J PC could be 0++ or 2++. 
As a new charmonium-like state, so far, it has been interpreted as P-wave charmonium state l60l l48l. 
the D*D* Q molecule [HI, HI], the cess state 

particular, Ref. [60| has presented the open-charm decay behaviors of the possible cc options, i.e. 
Xc,/(3P) (J = 0,1,2), using the quark pair creation (QPC) model. Only the decay properties of 
Xc2(3 3 P2) are consistent with the existing experimental data on the A"(4350). Besides, the hidden 
charm decay P(V(4350) — > J/ip(f>) [48[ via DD^ intermediate states favors the assignment being 
P-wave charmonium as well. 



631 and a mixture of scalar cc and D*D* [64| etc.. In 
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The mass of Xc 2 (3 3 P 2 ) in our calculation is 4383.4 MeV (NR) and 4352.4 MeV (MNR). Conse- 
quently, the quantum number of the V(4350) is reasonably favored as J PC = 2 ++ with the charmonium 
state x c2 (3 3 P 2 ) from the current exact mass prediction and strong decay behavior. To confirm this new 
member of charmonium family, more experimental data are required. We provide the corresponding 
radiative transition rates (see Table [6] for El channels and Table [7] for Ml) and the two-photo decay 
width P 77 (3 3 P 2 ) = 0.66 kcV for the further testing. 

4.6 Ml transition of J/tjj and ip' 

A well-known problem is evident in the decay rate of J/ip — > jrjc, which is that the predicted rate 
in the nonrelativistic potential model is about a factor of 2-3 larger than experiment in the previous 
calculations pll I29I] . Since this rate only involves the charm quark magnetic moment, this discrepancy 
leads to a surprise. When we adopt the wave functions generated from the full-potential Hamiltonian, 
this unexpected discrepancy seems can be reduced. 

The radiative transition rates of ip' from the previous perturbative calculation, both NR model 
and GI model, have showed serious disagreement with theory. Particularly, the predicted partial d ecay 
width of ip' — » jrjc was nearly one order of magnitude larger than the experimental observation [251 ] . 




The lattice QCD calculation [65( obtained -0' — > jt] c = 0.4 ± 0.8 keV, while its uncertainty was 
apparently large. Recently, Li and Zhao has interpreted this discrepancy of potential model via taking 
the intermediate meson loop (IML) into account (66|. The IML contributions bring down the GI 
amplitude to 2.051^75° keV, where the cut-off parameter played a crucial and sensitive role. Besides, 
the uncertainties from IML also due to the radiative couplings. We re-examined the NR model through 
the full-potential calculation, and found r(ip' — > 777c) = 3.1 keV (see Table [7]). In comparison with 
the perturbative results, the present approach may obtain a better wave function. For t/j' — >• jr]' c , our 
prediction is also in a reasonable agreement with the first measurement [33j of this decay channel by 
BES-III Collaboration very recently. 



By mean of a proper variational method, we calculate precisely the eigen value and eigen function of the 
Hamiltonian including the full interactions, which might maintain the subtleness of wave functions. We 
investigate the mass spectrum, electromagnetic transitions, leptonic decays and two-photon widths of 
the charmonium states. Comparing the results from the pure scalar and the scalar-vector mixing linear 
confinement, it is revealed that the mixture might be important for reproducing the mass spectrum 
and decay widths. Explicitly, about 22% vector component is obtained via optimizing the potential 
parameters. The scalar-vector mixing model requires less linear slope, and more in the spin-splitting 
than the pure scalar confinement. In our calculations, the Ml radiative transitions rates of J/ip and i[>' is 
reduced by the obtained wave functions actually, consisting with the experimental observations in order. 
The S-D mixing are explored in the higher charmonium states. To the newly found charmonium-like 
states, we perform a reasonable inspection of the possible cc assignments. The V(4160) and V(4350) 
are found in the mass spectrum as M(2 1 D 2 ) = 4164.9 MeV and M(3 3 P 2 ) = 4352.4 MeV, which 
strongly favor the J PC = 2 h ,2 ++ assignments respectively. With the available experimental data, 
the Z(3930) is favored as x c2 (2 3 P 2 ), while the V(3872) and V(3915) are disfavored by charmonium 
states. The corresponding radiative transitions, leptonic and two-photon decay widths of the V(3940), 
V(4160) and V(4350) have been calculated for the further experimental search and identifying the 
new members of charmonium family. 
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Appendix: Wave functions and matrix elements for Gaussian basis expansion 

Since the state wave function ipi m (r) are expanded in the Gaussian basis space, the Hamiltonian 
matrix elements and the radial wave function at the origin, as well as the overlap integral involved in 
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the transition matrix elements are explicitly analytic. We present these formulas used in the numerical 
calculation of the mass and decay width as follows. 
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